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J> i Abstract 

^ We study the Schrodinger equation which comes from the paraxial approximation of the 

Helmholtz equation in the case where the direction of propagation is tilted with respect to the 
i boundary of the domain. This model has been proposed in [12]. Our primary interest here is in 

the boundary conditions successively in a half-plane, then in a quadrant of R^. The half-plane 
p I ■ problem has been used in [llj to build a numerical method, which has been introduced in the 

^ ; HERA plateform of CEA. 

'. 1 Introduction 

The simulation of a laser wave according to the paraxial approximation of the Maxwell equation 
has been intensively studied for a long time when the simulation domain is rectangular and the 
direction of propagation is parallel to one of the principal axis of simulation domain(see for instance 
^16j and references therein). 

We are concerned here with the case where the direction of propagation is not parallel to a 
■ principal axis of the simulation domain, and cannot even be considered as having a small incidence 
angle with it. It may be crucial for example if we want to simulate several beams with different 
directions, possibly crossing each other in the same domain. This tilted frame model has been 
considered some years ago by physicists for dealing with beam crossing problems (see pi]), and is 
of particular interest in the framework of CEA's Laser Megajoule experiment (see [23| [25]). 

In [T2j [TT], a new model was proposed to deal with this case, and a numerical scheme was 
^ ' introduced and coupled with a time-dependent interaction model (this scheme was then used in the 
HERA plateform of CEA). We focus here on the theoretical study of this new model, of "advection- 
Schrodinger" type. 

This model is derived from the paraxial approximation, intensively used in optic models and 
in a lot of models related to laser-plasma interaction in Inertial Confinement Fusion (cf [HI HEi [13 
[20l[3l[21] and their references). According to laws of optics, the laser electromagnetic field may be 
modeled by the following frequency Maxwell equation (the Helmholtz problem): 

e^Aip + ip + ili^tetp = 0, Ic G V (1) 

where is the wave number of the laser wave in a medium corresponding to the mean value 
of the refractive index and vt = ^ + i^J' a complex coefficient : its real part u corresponds to a 



00 

o 

> 



*INRIA Rocquencourt, projet BANG, Domaine de Voluceau, BP 105, F 78153 Rocquencourt, Prance; email: 
marie.doumic@inria.fr 



1 



conveniently scaled absorption coefficient and its imaginary part n to the variation of the refraction 
index. 

The first basic assumption is that the wave length 27re is small compared to the size of the 
simulation domain (indeed, it is in the order of l//m and the simulation domain in in the order of 
some mm for the Inertial Confinement plasmas). Assume also that the light propagates according 
a fixed direction defined by the unit vector k . Let us denote the longitudinal and the transverse 
coordinates by z = Ic ■ k and X = Ic — {Ic ■ k ) k ; the gradient with respect to X is denoted 
V_L = V — k {k .V). Now, if we replace ■0 by u exp( ^ ), Equation ((T|) may be written as: 

e{2iutu + 2i k ■ Vu) + e^A±u = —e^-—ru. 

where Aj^ is the Laplacian operator in the transverse variable. Assuming that u is slowy varying 
with respect to the longitudinal variable, we can neglect the right-hand side of the previous equation 
(we will check this a posteriori ; see also [10] for more details). Therefore u satisfies the classical 
paraxial equation for wave propagation, which is a linear Schrodinger equation: 

i'k ■Vu + ^A_i_u + ii/tu = 0, 'x e V. (2) 

If the laser beam enters the simulation domain with a very small incidence angle that is, say, if the 
vector k is normal to a boundary of the simulation domain, the boundary condition on the entrance 
boundary {x = 0} is straightforward: it is simply the value of the solution, which is a data. On the 
face {y = 0} of the domain parallel to the vector k , to deal correctly with the boundary conditions 
one has to introduce a fractional derivative according to the x variable - see [HI HI El El US] for a 
mathematical approach and [18] for a physical approach ; see also Appendix lA.2[ or [2] for numerical 
treatment. _^ 

But, when the incidence angle between k and the normal vector to the entrance boundary ex 
is not small any more, what is of particular interest for the treatment of beam-crossing for instance, 
the classical approach is no more valid. 

The aim of this article is to analyse Equation jj]) with an arbitrary angle between k and ej, first 
in a half-space, and then in a quadrant, and to address the issue to find proper boundary conditions 
on the edges of the domain. The half-space problem has been studied and then used numerically 
in [11]; see also [7] for results of numerical simulations. To solve it, we first have to formulate the 
right entrance boundary condition and to show that the corresponding problem is well-posed. For 
this purpose, we consider a 2D problem but most of the ideas developed here may be extended to 
the 3D case. For the statement of the entrance boundary condition, one assumes that a fixed plane 

wave V'*'^ = tt*" exp(^-^-^) is coming into the domain where m*" is a given function of the variable 
y. Remark that for the Helmholtz problem, the boundary condition is classical and may read as: 

(e^+i'k ■li){iP-u'''e^^/') = 0, Ic e dV. 
o n 

where is the outwards unit normal to the domain. By using the definition 

. . 1 k . X . 

ip = uexp[ j 

e 

in this expression, the corresponding entrance boundary condition for Equation ^ may be written 
in a natural way as: 

{eTt ■V± + 2i'k7i){u-u''') = 0, Ic £ dV. (3) 
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In Section [21 for the sake of completeness and clarity, we recall (see [TT], Section 2) and detail how 
to analyse Problem ([2]) (El) in the half-space 



V = {{x,y) s.t. X > 0} 



which is the simplest case. First we recall a classical energy estimate in the space 
which implies uniqueness of the solution for any ut satisfying 



Then, in the case where is a stritly positive constant, we recall in Theorem [T] the exact analytical 
formulation of the solution, which is central in [TT] to build a numerical scheme, based on a splitting 
technique and on Fast Fourier Transform. We also develop here some corollaries on the regularity 
of the solution. 

We then have the tools to focus, in Section [3j on Problem ([2]) ([3]) in the quadrant 



which will make large uses of the solution on the half-space. 

The most delicate point is to find a proper condition on the boundary {y = 0}, which would be 
such that the restricted problem: 

1. is well-posed, 

2. admits for a solution the restriction to the quadrant of the solution in the half-space. 

Such a condition is called a transparent boundary condition. If the solution of the quadrant is not 
the exact restriction of the solution in the half plane, but is close to it in some sense, the boundary 
condition is called an absorbing boundary condition. 

For reasons which will appear in the sequel, the proper choice for such a condition is given by 
Equation (fTT]) . which writes: 



For a definition of fractional derivatives, see Appendix IA.2I or [22], ch. VI. 5. The derivation of 
this formula is similar to the derivation of the transparent boundary condition given for instance in 
[HI m El H] for the Schrodinger equation, using Fourier-Laplace transform: here, the operator k ■ V 
plays the role of time, and if we suppose ky = we would recover a "time" fractional derivative as 
in the previous references. 

Theorem [2] then reformulates Problem ((2l) (|3]l (fTTl) for functions whose derivatives may not be 
defined at the boundary. The main result of this paper is then given by Theorem [3l it states 
well-posedness and an explicit formulation for the solution to Problem ([2l) (|3]1 (fTTl) ^^^d analyses its 
relation with the half-space solution, showing that if the ray enters in the quadrant (i.e if /c -"n < 
) Condition (|17p is a transparent one, and only an absorbing one if it goes out of it {i.e if k - Ti < 0). 



inf-^i/( af) > 0. 



(4) 



V = {{x,y) s.t. x > 0,y > 0} 



dyU\y=o - A+{dx)my=o) = Vx > 



with a pseudo-differential operator defined by 
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2 Recall on the Half- Space Problem 



Here we address a problem where the simulation domain is the half-space 

p = {^.^>0}. 

Set 1e = {x,y) and define k = {kx,ky) to be a unitary vector of the plane which indicates the 
direction of the laser beam (we must of course have fc^; > 0). Assuming that holds and that // 
is a bounded function, we rewrite Problem dS])© under the form 

i k ■ Vu + -A_|_ti — fiu + ivu = 0, (5) 



(ieL>-2A:.7f)(M-M^")U=o = 0, (6) 

where we have denoted D = li ■ V± = ky{kxdy — kydx)- 

We first recall (See [TT], Proposition 1.) that if u G i7^(lR_|_ x R) is a solution to Problem ([5]) 
^ in the half-plane {x > 0}, we have the two following equivalent identities : 



J J 7 2 2\k ■n\ J V 2 2\k 



2i/|np + J \'k ■ 'n\\u\'^ = -Im{ J u{eD + 2ik ■ 7?)n^"). 



V To To 

The first estimate can be interpreted as the conservation of energy, the right-hand side being the 
incoming and the left-hand side the outcoming and the absorbed energy. >Prom the second estimate 
we deduce the following result. 

Proposition 1 Let {ieD - 2k ■ 'n)u^"- G L^(E). If u e H^(R^ x E) is a solution to the problem 
(E^ in the half-plane, it is unique. Moreover, we have the following stability estimate, with a 
constant C independent of u, fi: 



2u\u\'^ + J I'k ■Ti\\u\'^ <C J |(ieZ)-2T •T?)^™!^ 



V To To 

Remark: the stability result cannot be used to let the absorption tend to zero, since the norm on 
the left-hand side depends on the coefficient i'. 
In the sequel of this article, we assume that 

/i = 0, and u is a strictly positive constant (7) 

and we write 

'^kxg = icky{kxdy kydx)u -\- 2kxU 
Problem ©(IHl) may be written under the form: 

i{kxdx + kydy)u + ^{k^dyy — 2kxkydxy + kydxx)u + ivu = on x R, (8) 
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iekyik^dy - kydx)u\x=Q + 2kxU\x=Q = g. (9) 

We denote the Fourier variables of the variables x and y respectively by ^ and tj and the Fourier 
transform in x and y by J-x and J-y. Using Fourier transformation, the following theorem gives an 
explicit solution to this problem. Denote 



Ivy tl^y y 

Here and in the sequel, ^ denotes the principal determination of the square root (its real part 
is positive); it is crucial that v is strictly positive in order to define it without ambiguity. 

Theorem 1 Let S'{K) he the space of tempered distributions and g G S' . There exists a unique 
tempered distribution u G C^(1R+, 5^(R)) to Problem which is: 



Tyiu;x,rj) = 2.F,(g;7?) ^^^^ 

1 _L . / 1 _ 9£V? _L 9,-,.!^ 

It satisfies also: 



l + A/l-2^ + 2ii. ^ 



dx - R-{irf)jTy{u;x,ri) = 0. 

>From this theorem, we can also deduce the following corollaries. 
Corollary 1 If g £ H~^(M) then the solution u to Problem is in 

C^(E+,L2(ir))^ and the 

following stability inequality stands for some constant C not depending on the absorption factor v: 

On a general manner, ifge H'{R), seR, then u G C^(1R+, Hy^{R)) and we have the following 
inequality: 

\\u\\ , I 1 < Q rfsfTOV 

The space chosen for the solution is the most convenient because the stability obtained does not 
depend on the absorption coefficient z^. However, we can find existence and uniqueness of a solution 
in other spaces, provided other conditions on the initial data g. For instance: 

Corollary 2 1. If f"" ^^'''^ G LlCR) then u G L'^CR+ x R). 

2. If e L^i^) with s > then u G L^{R, Ll{R+)). In a general manner, if 

^"^'^vlT^S^t"^/'^'" ^ then u G {R, {R^)) . 

Notice that when i/ — > 0, TZe{R-.{iri)^=Q) values zero in a half-line, then with these assumptions u 
is not uniformly bounded according to the absorption coefficient u. 
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Proof. The proof is based on Fourier transforms of the problem, taken successively in y and in 
X. The factorization of the equation is the essential step of the proof, and will also provide us tools 
to define a proper transparent boundary condition for the quadrant problem. 

Let n be a solution of Problem ([8|)((9|) and v the extension of u by zero in the whole space: 
v{x,y) = u{x,y)lx>o- By introducing formally the function v in Equation ^ we get: 

ik ■Vv + ^Aj_v + iuv = (^{ik^ - ^{Ik^dy - kyd:r))u{0,y)^Sx=o + ^^(0, y)(5^=o- 
The meaning of the term dxu{0, y) is given by the entrance boundary condition ([9]) : 

ik • + |A_Lf + ivv = ikxg{y)5x=Q - ^ {kxdyu{0, y)6x=o - kyu{0, y)S'x=o) ■ 

Assuming that u £ C(1R+, 5'(]R)), we are allowed to take the Fourier transform of this expression. Let 
us define Pu{X,Y) as the polynomial which characterizes the differential operator of the equation, 
i.e: 



Pu{dxi dy) — i{kxdx + kydy) + ^{kydxx "^kxkydxy + k^dyy) + 



If we write uo{y) = u{0,y), the Fourier transform in y of the equation in v reads: 

Pu{dx,irj)^yiv;x,r]) = "^j (^^-^2/(5; ^) " ^p^-^j/K; f/) j '5x=o + -^^K; f?)'5Uo| • 



We can write 



ek 



PAdx,irj) = ^[dx- R+iiT]) ][dx- R-iiv) ) (11) 
where we define: 



l^y y ^ ' 

Thus: 

(^dx - R+iirf)] (dx - R-iirj) ]J^y{v;x,r]) 

^y ' ' 



kx \ ' 

^yia; n) - i-j^^^y^'^^'^ ^^=^ + ^yi^o; ri)^x=o- (12) 

We have reduced the problem in finding a unique acceptable solution for this ordinary differential 
equation (fT2]) . Its Fourier transform in x reads: 



(^i^ - R+iiT])^ (^i- R-{i'q)^J'xJ'y{v]i,r]) = '^^^y{9;v) - iijr^ - O^yiuo;v)- 

Since TZe{i£, — R±{i'r])) / 0, we can divide each side of Equation (fT2]) by -^Pu (because is an 
infinitely derivable function with at most polynomial growth at infinity, we can multiply it with a 
tempered distribution) and write: 

Vy{v;i,r]) = + 



— R+{i'rf) i^ — R-{i7f)' 
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where: 



"1/2 R^ ^irj) _ (i^) ' eA:2 (^,7) - R_ (iri) ' 

If G C\E,, we can verify immediately that: 

1 _ [ ^x.(l..>oe'";0 if7^e(0)<O 
^e-^ I -•^x(la=<oe'";e) if^e(0)>O 



(13) 



Applied to 9 = R+{iri) and 6* = R^{iri), since TZe{R+) = —lZe{R-.) > 0, we get: 

which proves that v G C(lR+, 5'(1R)) . Since v{x < 0, y) = 0, this equality implies a^^gl^) = 0- We 

take the limit for x ^ 0, and find !Fy{u] 0, = a^i^it])- It gives Formula ([TO]) and ends the proof of 
Theorem [D 



Notice that with Formula (fTOl) . we easily calculate the value of the derivative k ■ Vu. As soon 
as u is regular enough to have its Fourier transform in y decreasing rapidly, we can develop it 
asymptotically in e and z^, and find: k ■ Vu = 0{e + z^). It is a way to check a posteriori the validity 
of our first approximation. 

Proof of the Corollaries. 

To show the stability result of Corollary [H we use the following lemma: 

Lemma 1 There exists a constant C > 0, depending only on the geometry of the problem (\.e. of 
e, kx and ky) and not depending on v, such that: 

1 + 



\l + \ll-'-^V + 2ieS\^ 



Proof. Let us denote X = 1- '^r] and N = 2ei/§. Then: 

B{7]) = |1 + VX + iNl"^ = |1 + e'^VN-iXl"^ 

= 1 + + iV2 + 2{X' + N')hos{^ - ^I^). 
The cosine is nonnegative, so we get: 

B{rif > (1 + + iV2)2 > 1 + 

Since 1 + = 1 + a^(l - X)"^ with a = we study the polynomial P{X) = c(l + X^) - (1 + 
a2(l - where c = ^. We look for a c> such that V X e R, P(X) > 0. 

P = X^{c - a^) + 2a^X - (1 - c + a^), 
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and we pick c = 2a + 1 (for instance). We can take any C > 2y 1 + Tj^- ^ 
We now integrate in r] the square modulus of Equation (fTOl) . 
Since \e^-('v)^\ = e^<=(-R-(^'?))^ < 1, we get: 



CO 



We have proved the stabihty inequahty, and the fact that the solution u is in the space L'^ (R+ , L'^ (R)) . 
The constant C is of order The proof is the same for g G H^(Ji). 

In order to prove Corollary [21 we first show that if s > and ■^y^3,v)i^+\v\ )^ ^ L'iCE) then 

y/\ne{R-{ir}))\ I 

u G L^(1R,L2(]R+)). We have: 



We can apply Cauchy-Schwartz inequality: 

OO OO 

|u(x,y)|2dx<C2 / dxf /"dr?|.F,(5;r?)|2(l + |r?|2)V^4^-(^''))-V/"dr?(l + |r?| 



it gives, for some constant C > 0: 

\:Fy{g-^)\\^ + VW\?' 



This shows that u G L^(R, ^^(11+)) if we suppose for instance ^"^^^^^j^ [^^^^ ^ L'^r^i^). To gener- 
alize at i7™(R+), we replace J^y{g) by R-{i7])'^J^y{g;r]). 



3 Case of the Quadrant 

We still consider Problem dS])® but restricted to the quadrant {x,y > 0}. To find a transparent 
or an absorbing boundary condition (cf. Section [T|) on the boundary {y = 0}, we formally factorize 
the differential operator of the advection-Schrodinger equation as follows: 

= • V + |a^u + ii, = o = e'^{dy- A+id^)) {dy - A_ (a,)) . (14) 

where 



It means that Aj^{X)^ A^[X) are the roots of the polynomial Py{X,Y) considered as a polynomial 
of the variable X (see Appendix IA.2I for a recall on the definition of a fractional derivative, or ^22j . 
chapter VI. 5 for a complete construction). 
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The aim of this section is to show that the following condition: 

{dy - A+{d^,))u\y=o = V x>0 

is a transparent boundary condition if A;^ > and an absorbing boundary condition if ky < 0. The 
proofs will make appear more clearly why we were lead to such a choice for a boundary condition. 

Let us consider the following problem: 



i{kccd^ + kydy)U + -{kld'^y - 2k^kydly + kpl^)U + ii^C/ = V x > 0, y > 0, (15) 



ieky{k^dy - kydx)U\^=o + 2k^U\^^Q = 2kxg+ V y > 0, (16) 

dyU^y^o - A+{d,){Uiy=o) = Vx > 0, (17) 

where we take 5+ G i7^2(]R+), Supp{g+) C M*^. We can rewrite Condition (fTTl) under the following 
form: 



2iek k'^ 

iekx{kydx - k^dy)U\y=Q + ky{l + W 1 + — r2^t^x + 2ieu-^)U\y=o = V x > 0. 

\/ f^y rZy 



The first step is to give an extended meaning of the boundary condition JIZJ, for functions whose 
derivatives at the boundary {y = 0} may not be defined. The second step is to show that Condi- 
tion (fTTll . considered in this extended meaning, is either a transparent or an absorbing boundary 
condition for our problem, according to the sign of ky. 

Let us denote u the solution of the half-space problem jS]) ([9]) with an entrance boundary condi- 
tion equal to g = (7+ly>o: g S -?/~2(lR). It is given by Formula (fTOl) : 

Ty{u-x,rj) = _Ji(£i^)_^ei?.-(..)-. 



l + ^l-^rJ + 2^eu^^ 
We show the two following theorems. 

Theorem 2 Let g+ be a given function such that Supp{gj^) C and gj^ G i?~2(]R,_,_). Let U G 

-+s -- 

C^(lR+,ffy (1^+)); with s > 0, a function satisfying Equations / fT^jJ / fTfi]) withdyU{x,0) G Hx ^ 
and U{x,0) £ i?|(E+). Then U satisfies fT^) iff C/ satisfies: 

U{x,y)ly>o = :Fy'(^{k{r|)J^y{Uly>o■,0,v) + G{v)}e'^-^''^^^^^ (18) 

with: 



Rj^{irj) — R^{irj) ek^ Rj^{irj) — R^{iri) 

With this result, we are able to extend the meaning of Equation (fTTl) in the following way. 



Definition 1 We say that a function U G C^(1R+, Ly(lR+)) is a solution of Problem |73]) fT^) fT^) 
iff it is a solution of Problem [T5\} fT^) l[TS\} . 



9 



We are now able to state the main 

Theorem 3 Let g+ be a given function such that Supp{g+) C and g+ G H^2(K+). 
i) Problem ^l!^)1^10^)ilS\) admits a unique solution U G C^(IR+, Ly(IR+)). 

a) Letu G C^(]R.+ , Ly(R)) be the solution of the half-space problem f^^f^i with an incoming boundary 
condition g = g^\y>Q. We have an explicit formula for U in terms of u, with UQ{y) = u{0,y): 

U{x,y)ly>o = J'y^l^{K{v):^y{uoly>o;v) + G{v)}e'^~^'^^^^^^ (19) 



Where K = and G = 

Hi) If ky > 0, then U = M|y>o- 

iv) If ky < 0, and if we take g'^{y) = h{y — A) with A > 0, denoting u^,U^ the corresponding 
solutions respectively in the half-space and in the quadrant, we get: 

lim IKu"^ - ^^^)1?;>o||l-(r+,l2(R)) = 0. 
3.1 Fourier transform of the problem. 

Let [/ be a solution of Problem ([I5|)([l6l)([l7|, C/ G C^{R+, H^+%R+)) and V the extension of U by 
zero in the whole space: V{x,y) = U{x,y)lx>o'^y>o- We calculate what the problem means for V, 
and find: 

i'k ■ VV + + iuV = -ekr,kyU{0, 0)6r,=o5y=o+ 

\2 



ikx - ^-^{2kxdy - kydx))U{0,y)^ly>o6x=o + ^C/(0,y)lj,>oCo+ 

ek \ ek"^ i 

^^y - -^(.'^kydx - kxdy))U{x,0)jlx>oSy=o + -^U{x,0)lx>o6y^Q, 

which gives, with the two boundary conditions (fT^(fT71) : 



i'k ■ VV + |a_lF + ii/V = -ekr,kyU{0, 0)4=o'^y=o+ 



(^ikxgSx=o - ^{kxdyU\x=o5x=Q - kyU\x=o^'x=o)^'^y>o+ (20) 
^(1 - ^1 + ^dx + 2ieu'^)U\y=^5y=o - '-^{kydxU\y=o6y=o - kxU\y=o6'y^Q)^lx>o- 



To take the Fourier transform of this expression, we have to proceed carefully with the derivatives 
of U. Indeed, we can write: 

J^x{idxU)lx>o;C,0) =i^J^x{Ulx>o;C,0) -uio,o), 

and the equivalent formula for J^y(^{dyU)ly>o;0,r]) . On the contrary, one directly gets: 



/ / 2iek k'^ \ / 2ek k'^ 

^ ^ 2.ei.^C/|,=o)lx>o;e,oj =Jl- -^i + 2^ev^rx{U\y=o'^x>o). 
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We now take the Fourier transform of Formula ((201) . The terms depending on U{0,0) cancel one 
another, and we furthermore notice that: 



\ y y 



Finally, using the polynomial Pu{X,Y), we get 



Dividing by written in one of the form (fT^ or (fTTl) . the equation in F reads: 

^x^y ^ ' ^ ' - _ (^^) + _ (i^) + iry - (iO ir/ - ^_ (i^ ' ^ ^ 

where and are given by 



r(e) = o. 

3.2 Proof of Theorem [S 

Let us take the inverse Fourier transform of Equation (f2T]) . as in Section[2]for the half-space problem: 

If we multiply each side by la;>oly>o and use the fact that f5~ = 0, we obtain 

C/(x,y)l,>olj,>o = ,^j7'(^a-(7?)e^-(*'')"^l,>olj;>o, 

and denoting 



Rj^{i'r]) — R-{iri) ek^ Rj^iirj) — R-{i-r]) ' 

we have a~ {rj) = K (j])!Fy(U\y>Q]Q,ri) + G{ri), and get Formula (fTSll . 
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It proves that if U e C^^{R+,H.^ ''(1^+)) is a solution of Problem ([l5|) (US]) ([I7| , then it verifies 

Equation Conversely, let U € C^^{'R+,H^ (1R+)) be a solution of System (HH) (US]) (HE]) , and 

dyU{x,0) G Hx ^ and U{x,0) G Hx . Assuming no particular condition on the boundary {y = 0}, 
we carry out the same computation as in Section ISTTl and find: 



i£, — Rj^{iri) — R-{iri) ir] — A^{i^) iry — 
with the same than defined in (1221) and with: 



J'xl^idyU - A+{dx)U)lx>o;^,0 

^ ^^^^ A+ii^)-A_{iO ^ • 

We take the inverse Fourier transform and multiply it by lx>o^y>o, to find: 

Since U verifies also Equation (fTSl) . which writes: 

U{x, y)lx>oly>o = [or (ji)e''' (^"^"^ lx>oly>o, 

we necessarily have: 



The function JJ being continuous in y, we can take the limit when y tends to zero and obtain 
Tx^{ji~ =0 Vx > (see [l3] and the theory of Hardy functions). 

Since /?~ is the Fourier transform of a distribution null in R_ (see Appendix lA.ll) . it implies 
that /3~ is the Fourier transform of a distribution whose support is included in {x = 0}. Since 
f3~ € L^(IR), it implies f3~ = 0, which is equivalent to the Fourier transform in x of the transparent 
boundary condition. We have proved Theorem El 

Equation (fTTl) can be defined only for functions in i/2+* in the variable y. Theorem [2] allows us 
now to extend it by Condition (fTSjl . with Definition [1] given at the beginning of Section [3l 

3.3 Proof of Theorem [S 

We recall Formula (fTSl): 



U{x,y)ly>o = J^y^i^{K{v)J'y{Uly>o;0,v)+G{T])}e'^-^'^>j^ 

If we set Uoiy) = U{0,y)ly>o, it gives: 

Uo{y) = Ty^ [K(Ti)Ty(JJ^- ri) + G{ri)^ l^>o. (23) 

For these equations to make sense, it is sufficient to have U G C^(1R_|_, L^(1R+)). We prove the 
following lemma, from which we will easily deduce Theorem [3l 



12 



Lemma 2 i) Equation I118\) admits at most one solution U E C^(IR+, L^(IR_|_)); Equation 
admits at most one solution Uq £ L^(IIl). 

ii) The solution u G C^(]R,+ , L^(Il)) to Problem ^(3) (with g = g+ly>o) satisfies Equation / fl3|) . 
Moreover, U = u^yyo satisfies Equations /il5\} /il6]} . and if ky > it also satisfies Equation ^8\) . 



Proof. 

i) Uniqueness of a solution of Equation ([23]l in L^(Il) implies that of a solution to Equation ([T 
in C^(lll+, L^(1R,_|_)). Let us take g = and suppose that a function Uq G L^(R) verifies: 

Uo{y) = (K{r])Ty{UQ- r?)^ \y>o. 

Let Vo(y) = !Fy^^K{i])J^y{UQ;ri)^{y). The function Vq verifies the following equation, because 
Uq = Voly>o: 



Vo{y)=T-'(^Kir])J'yiVoly>o;r^)^ = J K{y 



s)Vo(s)(is. 





We separate Vq into y+ = Vblj;>o and V- = Voly<o. Since Vq G -^^^(11), V± £ the functions 

V± := J^y{V±;ri) belong to the Hardy spaces Ti"^^ (see Appendix I A. II for a recall on Hardy spaces). 
The resulting equation writes: 

y+(r?)(l - K{rj)) = Ul+ ^ =y+(r?) = -V-{ri). 

^ V L. 2eky , „. feg / 

The idea is to find an Ti.'^~ function on the left-hand side, equal to a TC'^'^ function on the right-hand 
side: since Tl'^^ Pi = 0, it will imply that both sides are null. We use the fact that the function 
can be extended analytically by a uniformally bounded function, respectively: on 



1 



-I 2tky ky 

]R2- ^ _ ^ > 0} if ky > 0, and on 1R,^+ if ky < 0. 

1. ky > : 1 — K{r]) is analytic and uniformly bounded in R^^, so (1 — K)V+ G 7{'^~ , and is 
equal to {-V-) G n^+, so V± = 0, so Vb = 0. 

2. ky < 0. We write the previous equation in the following form: 

V+{r^)=-{l-k{v))~'V.irj). 

(1 — K{r]))^^ is analytic and uniformly bounded in 1R^+, so (1 — K)^^V- G Ti'^^, and is equal 
to {-V+) G n"^-, so Vo = 0. 

ii) Let u G C^(1R+, Ly(lR)) be the solution of Problem ([HI)® in the half-space, with an entrance 
boundary condition equal to g € i?^2(R). We have seen in Section [2] that u verifies: 



i{x, y)ly>o = Ty^ (a-/2(r?)e^-(^'')") l,>o, 



with a-^i^ij]) = K{ri)Ty{u] 0, ??) + G{r]). 

This is almost Equation ([TH]) : we have a~ = — KTy{uly<o;0,rj). 



Let us distinguish the two cases, according to the sign of k. 



y 
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1. ky > 0. We have seen in Section [2] that u verifies: 

2J^y{gly>o;r]) 



J='y{u;0,ri) 



1 + 1/1-^^ + 2^^4 



Since the function , ^ can be extended analytically by a uniformly bounded 

V 

function in 1R^~, if g S -L^(1R+) then Ty{u;0,r]) G "H^^, so u{0,y) = u(0,y)ly>o- This shows 
that in this case, a~ = so u verifies Equation (fTSl) . 

2. ky < 0. Since u(0, .) G -L^(IR), K{r])J^y{uly<,o] 0, tj) G "H^"*", the function u verifies the following 
equation: 



n(0, y)ly>o = Ty^ (a" (r/))ly>o + ^^7' (^(t?)^^ (uly<o; 0, r?)) 1 



Since Er(?7)^y(ulj;<o; 0, ?]) G W^^, it is the Fourier transform of a function null in 1R+, so the 
last member of the equation values zero. This shows that m(0, y) verifies Equation ([23]) . 



Let us deduce Assertions ij, ii) and in) of Theorem [3] from Lemma [21 First, Lemma [2] proves 
the uniqueness of a solution to Problem (fT3]) (fTBIl (fTTl) . In the case where ky > 0, we have nothing 
more to prove: the unique solution of Problem (fTSl) (fT6l) ( fTTll is the restriction to the quadrant of the 
solution of the half-space problem jS]) 



Let us now take ky < 0. IfU exists, it verifies Equation ([23]) and we have seen that it is also the 
case of u, so U{0,y) = u{0,y)\yyQ. Replacing C/(0,y)ly>o by n(0, y)lj;>o in Equation ([T8]l verified 
by U, this implies necessarily that: 



We now have to prove that the so-defined function [/ is a solution of Problem ([T5]l ([TBI) ( [T71) . First, 
it obviously verifies Equation ([T8]l . so according to Definition [1] we have only to check Equations 
([ED and ([I6|). Let us define: 

V := (e^-^^^>{k{r^):Fy{uoly>o; v) + G{rj))} 



We have U = Vly>o^ so we can check Equations ([TSl) ([TBI) on V as well. We notice that V = u — r, 
with: 

Ty{r;x,rj) = k{i^)J^y{u^ly<^-i^)e''''^''^^^. 



The function u verifies obviously Equations ([151) ([16]) , so we have only to check them on the remainder 
r. Taking its global Fourier transform and multiplying it by P^, we find: 
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Since J^y{uoly<o; rj) £ H}'^ and ky < 0, the right-hand side is the Fourier transform of a distribution 
null for y > (see Appendix lA.ll for a proof). Thus the remainder r satisfies the Schrodinger 
equation (fTSj) in the quadrant. Furthermore, at the boundary x = 0, we have: 

J^y{r;0,r]) = K{r])Ty{uoly<o;r]), 

which is also the Fourier transform of a distribution null for y > 0. It proves that r(0, y > 0) = 0, 
so r(0, •)|s/>0 verifies the boundary condition (fT6]l . 

Estimate of the difference between the solution in the quadrant and the solution in 
the half-space 

It only remains to prove Assertion iv) of Theorem [3l For the sake of notations simplicity, we 
skip here the indices and denote u, U instead of w^, U^. Let us assume that ky < and denote by 
u the solution of the half-space problem and by U that in the quadrant. According to Equation 
(fT9l) . we have: 



U{x, y) = Ty^ [ e^-(*'?)^[^(77).^,(uoly>o; + G{j])\ ) \y> 



>0! 



and furthermore u verifies the following equation: 



u{x,y)\y^^ = Ty^[e''~'^"i^'-{k{Ti)Ty{u^^^^ 
so by taking the difference, we obtain the relation : 

(n-^7)ly>o = .F-^(e«-(*'')"K(r?).Fj;Mj,<o;0,r?)^lj,>o. (24) 

Let us now assume that g(%)) = /i(y — vl), ^ > 0, G //^2(]R,) and Supp{h) C M*^. We can 
write: 

:Fyiuo;v) = pB^L^^^==:FyiH;^)e~^^^, 

1 + Jl-2eJ|r/ + 2ii/e|| 

with a function H G L^(R). We then have uq = H{y — A), so ||uo1s/<o||l2 = ||-ff|j/<-AllL2, tends to 
when A +00. We use Relation ([24ll to estimate the difference between the half-space and the 
quarter plane solutions, since TZe{R-{irj)) < 0: 

||(U- C/)(x,y)lj^>0||L^(R+,L2(]a)) < ||i^jrj^(nol;/<o)||L2(R) < C\\H\y<_A\\L2(^^y 



Conclusion and perspectives 

A mathematical analysis has lead to an analytical form of the solution of the tilted paraxial equation 
in the case where the refraction index and the absorption coefficients are constant. We have proposed 
a convenient transparent /absorbing boundary condition for the problem on a quadrant, shown the 
well-posedness of the so-defined problem, and estimated the difference between the solutions on the 
half-space and on the quadrant. 

However, this boundary condition is non-local, which implies much difficulty for its numerical 
treatment (see [H] for instance); up to now, the choice made in [TT] to deal numerically with this 
boundary condition was to add an artificial absorption coefficient, as popularized in [5]. 
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We have also noticed that Formula (fTO]l allows us to check a posteriori the validity of the 
formal asymptotic derivation of Equation ^ from Helmholtz equation ([T]). Complete estimates to 
justify it rigorously is a direction for future research (see also [10] for a rigorous justification of the 
time-dependent problem). 

Acknowledgment. The author expresses very grateful thanks to Prangois Golse and Remi 
Sentis for their precious help, guidance, ideas and corrections. 



A Appendix 

A.l Fourier transforms of functions supported in R_|_ 

Background on Hardy classes 

For more details, we refer to [T3] . 

Definition: let h : C. The function h is said to be a Hardy function if it verifies the two 

following properties: 

i) h is analytic in 1R^~ = {u = a + ib, b < 0} 

oo 

ii) sup I |/if,| I2 = sup f \h{a + ib)\'^ da < 00 . 

b<0 6<0-oo 

The so-defined space is called the Hardy space and written as TL^" . 

Theorem 4 Let h : —^Ch£ Ti.^^ iff there exists a function f G L^(R+) such that: 

00 

V(j G R^', h{io) = [ f{x)e''^''dx. 





Fourier transforms of functions that are identically zero in R 

According to Theorem IH the following kinds of tempered distribution / are Fourier transforms 
of tempered distributions / supported in R+: 

1. fGH'^-iR): see Theorem HI 

2. There exists a polynomial P{iS,) such that / = P(i^)^(^), with g the Fourier transform of a 
tempered distribution g, Supp{g) C R+: indeed, we then have / = !Fx{P{dx)g)- 

3. There exists a Hardy function (j) G TL^" and a holomorphic function bounded on an open 
neighborhood of R^", such that / = (p^ . 

Examples: 



(a) the function ^' is a rational function with no pole in R^ U {00}, i.e. 

^ A- 



1=1 



with Xm{ai) > 0. 
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(b) There exists a>0, /3gIR, s>0 such that 



{l + i{az + P)y 

It is sufficient to prove this for ^{z) = jpp[^- On the one hand, we have: 

l + izGR'T: = izG (-00,-1) 

l + iz£K^ = zei(l,+oo). 

This imphes that ^ is holomorphic on (z* = e^^"^^ where Log means the principal 
determination of the Logarithm). 

On the other hand, ^' is bounded on a neighborhood of : if Xm(z) < we have 

mfzM = - = - < 1 

' ^ ^' \l + iz\' ((1- Jm(z))2+7^e(z)2)^/2 - ' 

4. The two previous cases imply that if 5 G Tl^~ and / = F{i^)g, with F a product of functions 
of the previous types, then / is the Fourier transform of a distribution whose support is in 
R+. 

A. 2 Definition of the square root of a diff"erential operator 

We recall here briefly how fractional derivatives are built, and refer to [22J, chapter VI. 5. for more 
details. For a > 0, we define 

where we have denoted x+ = xla;>o. It is an homogeneous function of degree a, which generalizes 
As for the polynomial case, one has 

va>i, ^ = i;-i. 

ax 

Taking this property as a definition, it allows us to define Ya by analytic continuation for o G C. 
For the case a £ — W*, we notice that the fact that the pseudo-function Pf{x\) is not defined is 
compensated by the poles of the F function on —IN (see also [22], chapter II. 3.) In particular, one 
has 

Yo{x) = sgn+{x), F-i = 5o, Y^i = 5^'-^\ / G M. 

Since derivation corresponds to the convolution by 5'q = Y^2, we can define {-^)"', for functions 
supported in R+, by the convolution by Y-i-a. It interpolates between , which is the convolution 
by ^0 = ^-1) and convolution by 5q = y_2. For the case a = ^, we thus have: 

1/2 



{^) / = *'-3/.*/ = i:(>-./.*/). 



We can also have a formulation by the use of Fourier transforms: 

Hy-i/2m = - 
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We obtain this formula by r(i) = and by passing to the hmit when e — in 

/ ^dx = 2 / e-'y -'^y dy = ^ . 
Jo Vx Jo Ve + «^ 

Finally: 

These two formulations give the two following equivalent ways to define properly the square root of 
the differential operator, which is, in our case, —ky — 2iekxdx — 2ieuk^. 

1. Using the Fourier transformation, we have: 

^x{{\l-kl - 2i€kxdx - 2i€vklu)tx>o\i^ = 



1 /2 

2. Using the expression {£) f = y_3/2 * / = £{Y-i/2 ^ /), we get: 

\/—ky — 2iekxdx — 2i€vk'iiu) = ^ ^ ~—j= dx [ — . ds. 

V " VTT J vx — s 
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